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ABSTRACT 



A cotpairison of data from all previous energy loss experirt^ts per- 
formed at tlie Naval Postgraduate School with new theoretical predictions 
is presented in this thesis. In addition, the data have been extended 
to include beryllium. With this extension, experiments have new been 
conducted on materials ranging in atonic number frem 4 to 82. Agreement 
between experiment and tlieory is excellent for the most probable energy 
loss. However, theoretical values for the half-widths of the energy loss 
distributions generally are small corpared to experirrent for thicknesses 
less than 5 g/cm^ and large for thicknesses greater tlicin 5 g/cm^. These* 
experiments \;ere conducted in the energy range from 50 to 100 MeV, with 
absorber thidenesscs from 0.7 to 5.0 q/cm^. These effects v.’crc found to 



be independeait of atonic nurnber or incident energy. 
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I. INTRODUCTION 



Fast electrons \'^ich traverse a layer of material will exit v/ith 
less energy than they possessed on entry. Che energy loss distribution 
can be characterized by two parameters, the most probable energy loss 
(Op) and the half-v/idth (HtV) , vvhich is the rPull width at half maximum 
of the distribution curve. Itese quantities are the basis for ccffrpari- 
son bebl^7een experirrental and tlieoretical results. Losses can be 
attributed to ttx> major effects, ionization and radiation (bronsstrah- 
lung) . Ionization losses for fast electrons have been calculated by 
Landau [1], Williams [2], and many others. Early experiments resulted 
in good agreauent with theory for the energy loss, but these results 
gave a hall-width which 'v^as greater than theo::^' at energies 1 MeV. 

Blunck and Westphal [3] modified previous theoretical calculations by 
introducing radiation losses. This additional energy loss is super- 
inposed on the ionization loss and results in a larger half -width and 
enhancement of the tail of the distribution. Other losses have been 
discounted as unirrportant relative to ionization and radiation. All of 
the theoretical calculations are limited to "thin" absorbers, i.e., 
those in v^ch is small compared- to the incident electron energy. 

Several measuronents of energy loss for high energy electrons through 
various materials have been conducted at the Naval Postgraduate School. 
All were performed within an energy range of 50 to 100 MeV. These 
include the works of Bumiller, Buskirk, Dyer, and Miller on aluminum 
[4,5]; Goocin.n on aluminum and copper [6]; DeLeui.1, Raynis, et al on 
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alxiiiinura, coppcjr, and lead [7,8]; and ^5^sbl:xx)ker and Sandquist on tin, 
lead, and gadolinium [9]. With the extension to beryllium, data are 
new available on various iraterials with atomic nutribers (Z) betaveen 4 
and 82. 

Experimental consistency with past experiments was insxored by 
rerunning several thicknesses of lead absorbers. Numerous and exten- 
sive checks were made on the conputer program (Appendix D) written by 
Mosbrooker and Sandquist [9] to insure that it was consistent with 
Blunck and Westphal [3] theory as modified by the beam folding tech- 
nique (Appendix C) viiich accounts for the energy spread of the incident 
beam. 

Cemparison of all experimental data taken at the Naval Postgraduate 
School with the theoretical results of Blunck and Vfestphal [3] is pre- 
sented in this thesis. Tne data on beryllj.un are eiiLirely new; the 
older data have been presented before, but the theoretical values given 
here are new and represent calculations made mth the corputcr program 
of Appendix D. 
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II . THEORETICAL CONSIDERATIONS 



The Blunck and tfestphal theory of the distribution for the energy 
loss of a bean of ironoenergetic electrons in passing through a layer 
of absorbing material assumes that the energy loss Q is small corpared 
to the initial beam energy, E^. Let W(Q)dQ be the probability of energy 
loss between Q and Q + dQ, and X be that portion of the loss Q due to 
radiation. Hence, the ionization loss is Q - X. Considering these tivo 
loss processes, the probability of energy loss is 



Q 

W (Q) dQ = Wj (Q-X) Wg 00 dXdQ (1) 

o 

\diere and Wg are the energy loss distri};)utions for ionization and 
radiation respectively. 

The Landau equation [l], as modified by Blunck and Leisegang [lO], 
is used for the energy loss distribution due to ionization. The distri- 
bution is expressed as a function of Landau's dimensionless paraxieter 
X as follavvs; 



Wj(Q)dQ = $(X) dX E 

n 



C Y 
n'n 



/ 



b^ + 




exp 



(X — X ) ^ 
n 

b^ + Y ^ 
'n 



dX (2) 



vhere X = + In 

The terms used in equations (2) and (3) are defined as follows: 

^n' '*'n' ^n constants given by Blruick and Leisegang [lO] and 
are used to fit Landau's distribution to a sum of gaussian functions. 

R is the absorber thickness in cm. 



E. 

1 

aR 



- 1.116 



(3) 
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The quantity "a" is a function of the atonic number Z, the atonic 
weight A, and the density p of the absorber; and 3(=^/c) of the 
electrons ; 



0*154 tA\ 

a = — MeV/cm. (4) 

A 3^ 

The quantity b^ is a correction to Landau's theory given by Blundc 
and Leisegang [lO]: 



j^2 ~ ^ 2 ni m 



aR 



In 



ra 



2 E. 
1 



I (1 - 3^) 
m 



(5) 



where the sunrnation is over the ionization potentials of the atonic 



electrons, and is the number of electrons with ionization potential 



Q is the average energy loss due to ionization (no radiation) for 
electrons of incident energy E^, and is given by Sternheimer [11,12,13] 
as f ollov-TS ; 



A t 

Q = ^ 
32 



B+0.43 + lnE^- 3^~C - “ log^^QpAic) 



m 



MeV (6) 



vhere t is the thickness in g/cm^, and the constants A^, B, C, X^, a^, 
and are parameters of the absorb^ material. These parameters for 
tin, lead, copper, beryllium, aluminum, and various other materials are 
listed in reference [13]. The parameters for gadolinium are not listed 
but were interpolated by Mosbrooker and Sanquist [9] . 

For Wg, the energy loss distribution due to radiation, Blunck and 
Westphal [ 3 ] give: 



Wg (Q) dQ = BoR (C/E^) ^ 



(7) 



where 



a = 1.4 X 10~3 £|. 



4/3 In + 1/9 
7 . 1/3 



an 






(8) 



and B is a normalizing factor = — . 



The distribution of total energy loss according to Blunck and 
Vvfestphal is obtained by putting equations (7) and (2) into equation (1) 
and performing the required integration. The result is the energy loss 
distribution for a single electron of incident energy For cortparison 

of theoretical and experimental values, this distribution function was 
used, with corrections to account for the finite energy width of the 
incident electron beam. This treatment is described in Section IV and 



I^pendix C. 
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III. EXPERIMENTMi PPDCEDURE 



Beams of electrons ranging in energy frcxn 53 to 95 MsV were obtained 
by using the LINAC of the Naval Postgraduate School. Electrons were 
elastically scattered £can the accelerator beam through 90° by a thin 
(1 mil) aluminum foil. The electrons passed through the varying 
absorber thicknesses and were energy analyzed by a 120° magnetic spec- 
traneter described by Kenaston, Luke and Sones [14]. 

The genercil experimental .arrangement viss similar to that used by 
Masbrooker and Sandquist [9]. (See Figure 1) . The energy of the inci- 
dent electron beam was determined by vising a nuclear magnetic resonance 
probe to measmre the field of the deflection magnet of the LINAC. The 
scattered electxons are counted on a ten channel plastic sciiitillatox 
counting system that consists of ten front counters and a single backing 
counter vhich is operated in coincidence with the front counters. 




Figure 1. Experimental Arrangement 
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Beryllim absorbers of tiiicknesses ranging from 0.738 g/aa^ to 
5.908 g/on^ were mounted on a remotely controlled ladder that could be 
moved up and davn vertically in the electron beam. This allowed us to 
collect data on one absorber thickness, then move the ladder and run a 
different thickness without having to enter the LIN?!C end station or 
change the LIKAC beam in any way. The absorbers were placed at about 
3 on frcm the scattering foil. 

A lead absorber of thickness 2.825 g/on^ was used to cortpare our 
experinental data with the results of Mosbrooker and Sandquist [9] . 

The results were in agreement within experinental accuracy. 

Since the electron beam scattered into the absorbing material was 
not monoenerg€itic , the energy distribution was measvured before and 
after passing through the absorbers. The zero thickness absorber dis- 
tributions found by passing the beam only through the scattering foil 
had half widtlis ranging from 0.21 MeV to 0.46 MeV, or an energy reso- 
lution of about 0.45%. Following the measurerient of these "zero peaks", 
the varying aJjsorber thicknesses v;ere moved into place and an energy 
distribution was determined for each absorber thickness . 

Because of the limited quantity of beryllium available and space 
requirements in the target chamber, it was not possible to place all 
eight absorber thicknesses on one ladder. Therefore, each time the 
beam energy was changed or the ladder arrangarent was changed, a new 
zero peak was neasured. This accounted for the possibility that the 
character of the beam might change. 

The data represent the number of scattered electrons detected by 
the coincidence counting system at the exit of the spectrometer as a 
function of the energy selected by the spectrcrneter . 
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Normalization was control leid by integrating the current frcm a 
Secondary Bnission Monitor. This instrument was mounted both up and 
doi'ffistream of the targets. There proved to be no significant difference 
in the width of the zero peaks due to the position of the monitor. Each 
data set (ten points) represents a given integrated charge, i.e., a 
certain number of electrons passing through the scattering foil. 
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IV. TREATMEITT OF DATA 



Tte data for the experiment were taken with the ten channel counting 
systCT as described in Section III. To reduce the data, account must be 
taken of the characteristics of the counting system. These include the 
energy spread of the counters, which is about 0.3% between counters or 
3% for the entire system. Accordingly, the energy seen by the individual 
counters must be determined and also correction must be rtade for the 
differing efficiencies of the scintillators and associated electronics. 

A ccnputer program has been writte:i v>?hich determines the correct relative: 
counter energies and also corrects for coxanter efficiencies. Background 
corrections, which v/ere quite small, \'7ere determined by periodically 
removing the scatterer and absoxjjars frau tiie bearri. 

The reduced data were then plotted and from these plots the experi- 
mental most probable energy losses and distribution half-widths Vv’ere 
measured. These results and their uncertainties can be found in Table I 
in Appendix A. 

Since the scattered electron beam is not monoenergetic and the Blunck 
and Westphal theory is. based on a monoenergetic incident beam, to ccrapare 
theory and experiment the theoretical model must be modified. This modi- 
fication is accomplished by representing the incident distribution as a 
histogram, each bin of v;hich has a definite energy, and then unfolding 
this distribution into the theoretical calculations carried out by an 
IK4 360/67 ccrrputer. The unfolding technique is described in Appendix C 
and the computer program used is in Appendix D. 
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The ccmputer program described in Appendix D also calculates the 
theoretical half-widths and most probable energy losses. The program 
normalizes the theoretical curve to the experimental data and plots 
both the theoretical curve and experimental data. These plots for 
beryllium are shOAn in Figures [3] through L26] in Appendix B. The 
measurable parameters were ccarpared and the final results are tabulated 
in Appendix A. 

The ccftparison of all data fron all previous energy loss experi- 
it^its performed at the Naval Postgraduate Sdiool Llt'JAC \-jcls acccuplished 
by normalizing each set of experimental data to the Blunck and Westphal 
theory as calculated by the canputer program described in Appendix D. 
The comparison of the theory parameters so obtained to the experimental 
data can be found in Tables II through VI in Appendix A. 
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V. - RESULTS AND CONCLUSIONS 



The theoretical predictions and experimental results for the most 
probable energy losses and half -widths for beryllium can be found in 
Table I, Appendix A. Our theoretical predictions and the experimental 
results frcm previous theses ccirpleted at the Naval Postgraduate School 
for aluminum, copper, tin, gadolini'um and lead can be found in 
Tables II through VI respectively in Ajpendix A. The theoretical data 
shown are for the Blunck and Kestphal theory as iTodified by the beam 
folding technique. In all t^les BW refers to Blunck and Westphal theor^^. 
Op represents the most probable energy loss, and is the distribution 
half -width. The column headed % represents the percentage by vdiich theo::^ 
exceeds experiment. 

Agreorent between the BluncJc and Westphal theory end our experi- 
mental results is excellent for the most pi'obable energy loss. The 
average difference between experiment and theory is less than 0.5%. 
However, agreament for half-width is good only for thicknesses less than 
1 g/cm^ and is fair for thicknesses between 4 to 6 g/cm^. This result 
in half-width ccmparison is in contrast to previous results where it 
was concluded that the Blunck and Westphal theory predicts satisfactorily 
for thin targets v^iere "thin" v;as defined as t (z) ^/^<13.5, viiere t is 
the absorber thickness (g/cm^) and z is the absoi±)er atonic number. 

The observed trend for HW (experimental) divided by HW (theoretical) 
as a function of t (absorber thickness) is plotted in Figure 2. 
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FIGURE 2. Trend Plot of HW /RW.^ 

exp' th 



4 5 

vs Thickness 



t (g/cm^) 



A typiCcil ccrnpcrrson hct',.'CGn cr^poriincntd end thscrsticsl cur'/es for 



HW plotted as a function of t (absorber thickness) is presented in Figure 3. 




Tteoretical Half-widths vs Thickness 
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As can be observed from these plots the theoretical half -width 
tends to be srraller than experiment until one gets to thick absorbers 
(approximately 5 g/cm^) at 'vdiich point it gets larger than experintant. 

We have considered two possible effects vdiich might explain 
this crossing of experiment and theory takes place. One effect is that 
this trend is caused by an assunption that is made in the theory. This 
assumption is that of constant incident energy (E^) , i.e., that the 
energy loss is small cotipared to the incident energy of the electrons. 
However, in tlrlck absorbers, this assumption is not valid, therefore 
the theory tends to overestimate the energy loss due to radiation 
(bronsstrahluiig) . TlrLs overestimate will cause an increase in the dis- 
tribution tail and hence increase the distribution half -width. This 
trend in thick absorbers is not surprising in that the theory was 
originally designed for thin absorbers. The fact that the agreement .is 
excellent for wMle not for half-^i^idth is also not surprising since 
an overestimate of this type would not affect the most probable energy 
loss significantly. Although this theoretical inaccurac^^ is undoubtedly 
a cause of some discrepancy in agreement with experiment, it is doubtful 
that it is the cause of all the disagreement. 

The second possibility which we have considered affects the experi- 
mental results. Due to multiple scattering in thick absorbers the possi- 
bility existed that the size of the spectraneter opening v;as excluding 
the lew energy portion of the distribution. However, calculations of 
the projected mean square sc:attering angle as outlined in Jackson [15] 
shewed this angle to be approximately 0.30° while the spectraneter 
opening v;as approximately 0.90°. Therefore, it v;as determined that 
multiple scattering was not a factor in the disagreement of theory amd 
experiment. 
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Discxepancdes in the theoretical predictions frcm oxir program and 
previous results of Goodwin [6] and DeLeuil and Raynis [8] were noted. 
An atterpt was made to resolve these differences. The parameters from 
Blunck and Westphal [3] were inserted into our ccrtputer program and 
without the unfolding technique described in Appendix C v.e were able to 
duplicate closely the Blmck and Westphal distributions. On this basis, 
it is assuned that our theoretical results are the rrore accurate. 

Contrary to previous theses, the thickness at which theory and 
experiment began to exhibit a noticeable discrepancy was found not to 
be dependent on the atonic nmiber of the material. 

It should be mentioned that no atteirpt was made in this experiment 
to observe the t'ery low energy tail of ber^’^llium (approximately at 15 
MeV) where in previous experiments secondary effects had caused a peak 
in the absorber distribution. 

In light of what has been found in this thesis, further work is 
justified in several areas. First of all, the lo// energy tails of the 
various absorber distributions should be investigated to find out if 
secondary effects do cause a low energy pea}c in sane distributions. 
Secondly, an effort should be rtade to improve the efficiency of the 
current carputer program. Finally, an attanpt should be made to extend 
the theory to thick absorbers vhere the energy loss is not small with 
respect to the incident energy, and also to include the effects of 
secondary processes. 
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APPENDIX A - TABLES 



Table I. Energy Loss Distribution Characteristics of Berylliun 



Be(z = 


4) 


Qp(MeV) 


IM(MeV) 


E^ (MeV) 


t(g/c3n2) 


BW 


Experimental 


Q. 

*o 


BW 


Experirrental 


% 


52.89 


0.742 


' 0.98 


0.98 


± .02 


-.30 


0.33 


0.36 


± .02 


-8.3 




1.479 


2.04 


2.00 


± .04 


+2.0 


0.53 


0.56 


± .04 


-5.3 




2.209 


3.05 


2.96 


± .04 


+3.0 


0.67 


0.77 


± .04 


-13.0 




2.961 


4.17 


4.09 


± .08 


+2.0 


0.90 


1.07 


± .07 


-15.9 




3.673 


5.29 


5.19 


± .05 


+2.0 


1.17 


1.25 


± .06 


-6.4 




4.415 


6.36 


6.24 


± .13 


+2.0 


1.40 


1.46 


± .09 


-4.1 




5.179 


7.46 


7.40 


± .08 


+ .80 


1.76 


1.76 


+ .10 


0 




5.908 


8.69 


8.38 


± .12 


+3.6 


2.09 


2.19 


± .14 


-4.6 


74.78 


0,738 


1.02 


1.01 


± .03 


-1.0 


0.46 


0.47 


± .04 


-2.1 




1.479 


2.04 


2.01 


± .04 


+1.5 


0.56 


0.61 


± .04 


-8.2 




O ono 


3.11 


O A /I 


j- AC 


+2.3 


(\ 1A 

» i '.X 


0.79 


± . 05 


-6.3 




2.941 


4.15 


4.15 


± .14 


0 


0.96 


1.15 


± .08 


-16.5 




3.673 


5.29 


5.15 


+ .06 


+2.7 


1.17 


1.39 


± .06 


-15.8 




4.435 


6.39 


6.28 


± .10 


+1.7 


1.46 


1.61 


± .10 


-9.3 




5.179 


7.46 


7.42 


± .13 


+ .50 


1.76 


1.83 


+ .08 


-3.5 




5.908 


8.68 


8.51 


± .09 


+2.0 


2.12 


2.24 


± .11 


-5.3 


94.64 


0.738 


1.02 


1.05 


± .04 


-2.7 


0.58 




± .05 


-1.7 




1.479 


2.04 


2.06 


± .04 


-1.0 


0.68 


0.66 


± .04 


+3.0 




2.209 


3.11 


3.09 


± .04 


+ .60 


0.83 


0.84 


± .06 


-1.2 




2.941 


4.15 


4.06 


± .06 


+2.7 


1.06 


1.22 


± .08 


-13.1 




3.673 


5.29 


5.19 


± .06 


+1 . 9 


1.23 


1.35 


± .08 


-8.9 




4.435 


6.39 


6.30 


± .12 


+1.4 


1.53 


1.58 


± .12 


-3.2 




5.179 


7.61 


7.47 


± .12 


+1.9 


1.82 


1.93 


± .10 


-5.7 




5.908 


8.69 


8.63 


± .14 


+ .60 


2.16 


2.29 


± .11 


-5.7 
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Table II. Energy Loss Distribution Characteristics of Aluminum 



A1(Z = 


= 14) 


Qp(MeV) 


HW(MeV) 


(MeV) 


t (g/on^) 


BW 


Experimental 


% 


KV 


Experimental 


% 


53.57 


0.730 


1.05 


1.03 


+1.9 


0.44 


0.49 


-10.0 




1.441 


2.12 


.2.16 


-1.9 


0.62 


0.75 


-17. 0| 




2.146 


3.22 


3.18 


+1.2 


0.88 


1.02 


-14.01 




2.859 


4.29 


4.21 


+1.9 


1.20 


1.32 


-9.1 


53.44 


5.574 


8.86 


9.24 ± .13 


-4.1 


3.12 


4.60 ± .30 


-32.0 

1 


74.63 


0.730 


1.05 


1.04 


+1.0 


0.54 


0.51 


+5.9 




1.441 


2.12 


2.18 


-2.7 


0.72 


0.91 


-21.0 




2.146 


3.22 


3.17 


+1.6 


0.95 


1.08 


-12.0 




2.859 


4.37 


4.26 


+2.6 


1.25 


1.38 


-9.5 


74.47 


5.574 


8.86 


9.33 ± .25 


-5.0 


3.07 


5.7 ± 1.2 


-45.0 


96.93 


0.730 


1.10 


1.17 


-6.0 


0.64 


0 . 66 


-3.0 




1.441 


2.16 


2.16 


- 


0.81 


1.01 


-20.0 




2.146 


3.22 


3 24 


-0.6 


1.03 


1.19 


-13.0 




2.859 


4.37 


4.32 


+1.2 


1.31 


1.48 


-11.0 
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Table III. Energy Loss Distribution Characteristics of Copper 



Cu(Z = 


= 29) 


Qp(MeV) 


HW(MeV) 


E (MeV) 


t (g/on^) 


m 


Experimental 


% 


H4 


Experimental 


% 


52.84 


0.711 


0.92 


0.94 


-2.1 


0.46 


0.52 


-12.0 




1.423 


1.89 


1.87 


+1.1 


0.69 


0.84 


-18.0 




2.134 


2.90 


2.87 


+1 . 0 


1.04 


1.16 


-10.0 




2.845 


3.87 


3.80 


+1.8 


1.54 


1.75 


-12.0 


53.63 


4.295 


6.10 


6.23 + .13 


-2.1 


3.70 


4.25 ± 1.20 


-13.0 


74.76 


0.711 


0.92 


0.92 




0.56 


0.66 


-15.0 




1.423 


1.89 


1.91 


-1.0 


0.78 


1.07 


-27.0 




2.134 


2.90 


2.92 


-0.7 


1.12 


1.43 


-22.0 




2.845 


3.96 


3.98 


-0.5 


1.61 


1.98 


-19.0 


76.02 


4.295 


6.10 


6.28 i .13 


-2.9 


•3 TO 
1 ^ / O 


4.55 ± -50 


-17.0 




5.726 


8.65 


8.38 ± .25 


+3.2 


10.37 


10.1 ± 2.0 


+3 . 6 


94.30 


0.711 


0.91 


0.95 


-4.2 


0.67 


0.74 


-9.5 




1.423 


1.94 


1.94 


- 


0.87 


0.99 


-12.0 




2.134 


2.96 


2.93 


+1.0 


1.22 


1.47 


-17. C 




2.845 


3.99 


4.00 


-0.2 


1.71 


1.75 


-2.3 


96.66 


5.726 


8.82 


9.20 ± .50 


-4.1 


10.50 


16.4 ± 4.2 


-37.0 
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Table IV. Energy Loss Distribution Chciracteristics of Tin 



Sn(Z = 


= 50) 


Qp(MeV) 


HW(MeV) 


Ei(MeV) 


t(g/cni2) 


BW 


Experimental 


ct 

"o 


BW 


Experimental 


Q. 

'O 


52.53 


1.485 


1.83 


1.78 ± .05 


+2.8 


0.75 


0.78 ± 


.10 


-3.8 




2.970 


3.84 


3.86 ± .15 


-0.5 


3.03 


3.02 ± 


.40 


+0.3 




4.455 


6 . 56 


6.35 ± .30 


+3.3 


16.70 


11.30 ± 


2.0 


+48.0 


75.00 


1.485 


1.81 


1.77 ± .07 


+9.0 


0.73 


0.83 ± 


.10 


-12.0 




2.970 


3.93 


3.86 ± .15 


+1.8 


2.56 


2.96 ± 


.42 


-14.0 




4.455 


6.56 


6.65 ± .30 


-1.4 


16.69 


15.50 + 


2.0 


+7.7 


94.40 


1.485 


1.88 


1.87 ± .10 


+0.5 


0.83 


0.88 ± 


.15 


-5.7 




2.970 


3.93 


3.96 ± .20 


-0.8 


2.63 


2.72 ± 


.45 


-3.3 




4.455 


6.56 


7.00 ± 1.0 


-6.3 


16.75 


16.20 ± 


4.0 


+3.4 



Table V. Energy loss Distribution Characteristics of Gadolinium 



Gd(z = 


= 64) 


Q^(tfeV) 


HW(I4aV) 


Ej_ (MeV) 


t (g/an^) 


BW 


Experimental 


Q. 

■o 


BW 


Experimental 


% 


53.18 


0.814 


0.97 


0.92 + .08 


+5.4 


0.46 


0.52 ± .07 


-12.0 


52.80 


1.610 


1.94 


1.92 ± .10 


+1.0 


0.91 


1.00 ± .15 


-8.9 




3.221 


4.24 


4'.16 ± .30 


+1.9 


5.52 


5.75 + .72 


-4.0 


75.00 


0.814 


0.97 


0.97 ± .05 


— 


0.51 


0.59 ± .08 


-13.0 




1.610 


1.97 


1.92 ± '.10 


+2.6 


0.96 


1.00 ± .15 


-4.0 




3.221 


4.24 


4.09 ± .30 


+3.7 


5.55 


5.80 ± .80 


-4.2 


91.80 


0.814 


0.97 


0.97 ± .05 




0.60 


0.64 ± .08 


-6.2 




1.610 


1.97 


1.90 ± .10 


+3.7 


0.96 


1.10 ± .18 


-13.0 




3.221 


4.20 


4.13 ± .40 


+1.7 


5.56 


8.10 ± .95 


-31.0 
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Table VI. Energy Loss Distribution Characteristics of Lead 



Pb(Z 


= 82) 


Qp(MeV) 


HW(MeV) 


Ei (MeV) 


t(g/cm2) 


BW 


Experirtventail 


% 


BW 


Experimental % 


53.85 


0.706 


0.79 


0.81 ± .05 


-2.5 


0.54 


0.60 ± .06 -10.0 




1.412 


1.66 


1.45 ± .08 


+14.0 


0.97 


1.23 ± .12 -21.0 




2.118 


2.55 


2.64 ± .10 


-3.4 


2.11 


2.65 ± .25 -20.0 




2.825 


3.57 


3.66 ± .12* 


-2.5 


6.45 


6.25 ± .50* +3.2 


74.74 


0.706 


0.83 


0.82 ± .05 


+1.2 


0.78 


0.82 ± .05 -4.9 




1.412 


1.68 


1.80 ± .10 


-6.7 


1.23 


1.47 ± .07 -16.0 




2.118 


2.66 


2.59 ± .12 


+2.7 


2.37 


2.24 ± .19 +5.8 




2.825 


3.66 


3.70 ± .15* 


-1.1 


6.52 


6.25 ± .85* +4.3 


91.37 


0.706 


0.90 


0.87 ± .07 


+3.5 


0.93 


0.78 ± .08 +19.0 




1.412 


1.80 


1.77 ± .10 


+1.7 


1.37 


1.52 ± .19 -9.8 




2.118 


2.70 


2.75 ± .15 


-1.8 


2.53 


2.74 ± .64 -7.7 


— _ ■ . i 


2.825 


3.69 


3.84 ± .22 


-3.9 


6.95 


8.80 ± .70 -22.0 



*Experimsnt data average of DeLeuil/Raynis and Mosbrodker/Sandguist 
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Figure 11 
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Figure 13 
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Figure 14 
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Fig'jre 17 
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Figace 20 
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Figure 21 
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Figxire 24 
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Figure 26 
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Figure 27 
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Figure 28 
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Figore 29 



APPENDIX C - BEAM POIBING TEOMIQUE 



The theory of Blunck and Vfestphal [3] assunnes that the beam of 
electrons striking the absorber is roonoenergetic. The beam of elec- 
trons produced by the Naval Postgraduate School LINAC, or any other linear 
accelerator, for that matter, has a finite energy spread about a maximum 
or most probable energy point. The fact that monoenergetic electrons 
are not available to strike the absorber must be taken into account in 
ccitputing theoretical predictions if a meaningful coirpariscn with 
experimental results is to be made. This has been done here by a 
technique termed "beam folding" . 

Beam folding is acoorplished by a number of well defined steps. 

Energy distribution curves may be approximated by histograms. These 
histograms consist of a series of "bins" of area W (Q) AQ, vhere Q is 
the energy loss and W(Q)AQ is the probabiliby of loss between Q and 
Q + AQ. Thus, each bin has an "address", Q, on an energy coordinate 
scale. To acocnplish beam folding, the beam distribution must be known. 
This is observed experimentally and approximated in the corputer by a 
gaussian curve of appropriate half -width. The energy at vhich the 
maximum occurs is established by the energy of the beam, and the inagni- 
tude of energy loss incurred by elastic scattering of the beam as it 
impinges on the thin aluminum scattering foil prior to striking the 
absorber. The width AQ is then selected for the predicted distribution. 
This must be a small, but finite number vhere numerical techniques are 
to be used. This same width is used to break the beam distribution into 
a histogram. This is illustrated in Figiore 30 (a) . The reason the same 
width is used for both distributions is to facilitate ccnputer prograitrriing . 
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Figure 30. Beam Folding Distributions 
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Each bin of the beam distribution is now treated as a monoenergetic 
beam with energy coimensurate with its center location on the energy 
scale. Each bin also has a definite magnitude, or weight, with the magni- 
tude of the center bin being unity. The formulae of Blunck and t-festphal 
is now applied to each of these beams and an absorber distribution curve 
results for each, with a maximum arrplitude proportional to the height of 
the appropriate beam distribution histogram. Each of the absorber dis- 
tribution curves thus obtained may be thought of as being plotted and 
added to previously determined curves, using the energy of the electrons 
as a correlation point, as depicted in Figures 30(b), (c) and (d) . This 
is accomplished on the computer by adding the contents of each bin of the 
same address and plotting the cumulative total. The bin of address E is 
depicted in Figures 30(b) , (c) and (d) . Note that the beam distribution 
for histogram E4 is centered over the maximum point of the distribution. 
This is important as a false picture could easily be presented if the 
histogram were not symmetrical as the distribution would then appear 
skewed. 



54 



’>oo 



APPENDIX D - COylPUrER PPOGRAM 



C PLiRPDSci: This PPCG»TAM accepts parameters FP('M energy loss 

C EX?hRIME;NTS CEMOUCTEU UM THE 1. IHAG AMD ''.'ILL CCMPOTE 

C TMEG^ET ICAL HALFvUDTH AND MOST PRCJBAGLE LNEPGY LOSS AMD 

C VdLL PLOT ROTh THE EXPERIMENTAL AND THEORETICAL ENERGY 
C DISTRIBUT ION FDR THE EXPERIMENT. 

C 

C SCOPE: THE PROGRAM IS SET UP TO MAKE COMMUTATIONS RELATIVE 
C TO EXPER IME nITATI ON V,' I TH A LOMI NO M , COPP E R , T I N , G ADQL I N I UM 
C LEAL AND BFRYLLIU’'. 

C 

C MErnOO: THE PROGRAM PLOTS A THEORY CURVE BASED ON THE FIT 

C OF GAUSSIAN CURVES PERFOP.MED BY ul.UNCK A^jD NESIPHAL EUR 
C THE PkEUICTE'i ENFRGY DISTRIBUTION CURVES. THE IMTEGRA- 
C TIO'I OVE”- Ei'lIrRGV LC'SS iNwICAlED IN THE AFOR £ MF N T I ON ED 

C PAPER IS .^ekFCRMFO ON THE COMPlMER USING 3? POINT v,AlJSS 

C OUABRATURE. 1 NO GRAPHS \kV PLOTTED; ONc: REPRESENTS IHE 
C R-'n' PR-DICriON WITH A MON ENER Gtl I C t3EA> OF ELECTIONS AND 
C THE OTHER, WHICH ALSO HAS EXPERIMENTAL RESULTS THEREIN, 

C REPRESEiNTS THf S A'l E THEL.RY WITH A POL YE NE RGE T I C BEAM 
C 0 I S 1 R 1 BU T I ON FOL HE D IN. 

C 

C THIS PROGRAM IS WRITTEN USING FORTRAN IV LANGUAGE AND THE 
C NPS COMPUUR PaCILITY PLOniNG PACKAGE. Tl;TAL RUNNING 

c time is a euuCtion up thickness, atomic number, aNl 

c halfvidth. 
c 

implicit REAL«? { A-H,0-Z ) 

ESTABLISH TlTI.rS FOR PLOTS. 



P.FAL=:g3 TKA) 

CATA Tl/' 

REAL-B T3(A) 

DATA T3/ ' 

REAL^ti TA(A) 

DATA TA/'T: 
RcALYP T'3(A) 

DATA TB/' 

REAL'^.” T6(A) 

DATA To/ 'THEORY; 
RFAL-:=B T7(A) 

C A T A 17/' 

REA1>.^8 T3(A) 

CATA TJ/' 

Rt AL -H TG{ 4 ) 

DATA T-»/' 

REAL^^^B T 10(A) 
DATA rij/' 

REAL* a Til (A) 
CATA Til/' 
DI.-'ENSION l(pG), 
IWUFI 3c . ) ,XUF(3 j - 
EXTERNAL FCT,pFA 
DC 12.0 J = l,4 
El J)=.j 
123: Wr(J)= 

DC 12^1 J=1,2lj 
EXX{J)=l 
bXY( J)=N 
12 El IXY(J)=. 

0012. ■•2 J = 1,3'J0 
V!CT { .J ) =M 
Ha { J ) = 

X X { J ) = 

WUF { J ) = ) 

AD2 XUF( J ) = ) 



ENERGY D I ST PI DOT I ON 
ALUMINUM ABSiTRBER 
GM/CM LI: 

UNFOLDED THEORY 

EXPEP IMl-.NTAL data 
LEAD ABSORBER 
T-IN ABSORBtP 
GADOLINIUM ABSORBER 
COPPER -ABSORBER 



• / 

• / 

MEV / 

' / 
X* / 

' / 

• / 

• / 

' / 

/ 



BPRYLLIU-1 absorber 

NF (S'), V.MT { 3 t ) , HW ( B( R ) , X X { 3C ' ) , 

) ,PXX(2(' ') ,EXV(2- (') , IXY{2/u) , T"(A) 
M, '! I GN , l.'RAD 



IP 
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f 





onoooon ono non ooono oooo ooooo oooo 



c 


READ IN 


PARAMETERS FOR EXPERIMENT. 




c 


EO 


BE\M E MERCY 




c 


1 


ATCMIC number 




c 


A 


ATOMIC HEIGHT 




c 


T 


TARGET THICKNESS; GM/SQ CM 




C 


CZERU 


SPECTPOMEl FR INDICATION OF LOCATION 


OF BEAM 


c 

c 


CT SNJR 


DISTRI'fUTION (MEV) 
r, UMBER TO 'aHICH EXPEPiriMTAL COUNTS 


WILL BE 


c 

c 


HI 


NOKI'AL I ZED 

HALFWIDTH (F OEAM DISTRIBUTION 





RFA0(5, 11 ) ta, Z,A,T , Q7FRU , CTSIJUR , H T 
'WRIl^ rri3tZtA>rf l71'KOtCIl><l'RTHl" 

21 FCRMAl ( / ,2X , > !: : ',F?.2,‘ Z: ',FA.^, ' A: *,F5.3, 

1’ T: ',F 7 .a,' jZEKU: ',rF.3,' CTSNJk: ',F6.1, 

2 • HI : • , Ft>. 3,/ ) 

Divnu ZtRiJ PCAK I, '.TO "UINS" OF HIOTH DFLX. BIN hlDTH 
INCFEASr.S WllH TApOLF TH1CKNESS« 

0ELX=.i 

IF( T,LI . 1. ) DE LX = . •. 6^’^T 

CAVE IS THL RFCOIL LOSS OF FLECTknNS SIRIKlNG ALU'-'INUM 
NUCLEI 1 iF rut SCATTtklNf. FOIL. IONIZATION LOSS IN 
SCATTERFf’ IS NEGLECTED. 

C A V L = F B E B / ( 3 i . * 2 7 . ) 

CURB IS- IMF CORRECTION TO RE AITDED TO b XCFR I MENTAL ENERGY 
TO CORRECT FOf->, S EEC T RCME TL R ERROR. 

CORP. = QZ£RU-CAVF 
1 2 FORMA T ( F ( 1 X , F8 . 3 , I X , I 6 ) ) 

READ IN 10 2 E EXPERIMENTAL POINTS; EXX IS ENERGY (MEV) 
ANO IXY IS COOiNI S. TliERt. MUST FIE An UATa CakDS ARO 1 HE 
LAST FXX MUST BE ZERO. 

REAOIS, 12) (EXX( J) , IXY( J) , J = 1 ,2G0| 

establish parameters fop SCATTERING FOIE. 

FfJIL-.( }68A 
Z"aiL=l 3. 

KSIG=n 
00 1 J--=1,2.F 
I EXY( J)=IXY( J) 

P=) 

El IS VALUE OF BEAM ENERGY INCIDENT ON TARGET. 

El =EB-OAVE 
T E -IP^F 1 

THE FOLLOWING PORTION C'F THE PROGRAM UNECLOS BEAM UISTRI- 
iUJTlON. ICENT IS AN INTEGER VALUE ASSIGN'ED TF THE 
HISTOGKA*-' CENTER-D OVEP THE PEAK OF TrU: BEAM niSFRIPU- 
TION. EACH HISTOGRAM IS DENOTED BY ITS ENERGY E AND 
MAGNITUDE RELATlVt TO THE PEAK OF THE D I S T R 1 ^UT 1 ON , WF. 

ICENT= ( HT + or-LX ) /OFLX 
X^EI- ICENT^'-DFLX 
? X=X+DFLX 

IFIX.GT. (El+HT ) ) GOTO 1C03 
MOm+1 
F ( M ) = X 

WF( M) = F3GaM( E I , HT, X ) 

FS'GTO^.jM 

10^3 WRIT=f(a, 82) CAVE ,HT ,FOI L , ZFOI L , El 
O'’ 1 io L= I , ■; 

DX = n 

N=.: 
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K‘05 .\ = N+l 
C 

C USIKG EACH histogram AS A "BEAM", THE FOLLOWING PORTION OF 
C THE PROGR\M COMPOTES AN ENERGY D 1 S 1 K I BUT I ON FOR EACH 

C OF THESE ''BEAVS"' AND STORES THEM lO "BINS" LAbELEO BY 

C ENERGY, XX, THE CUMULATIVE AMOUNT OF ELECTRON CoUNTS IN 

C THE BIN IS WOT. Q IS THf ENERGY LOSS OVtER WHICH THE B-W 

C PRECICTIJG IS IMTEGRATEU AND OX IS TrIE AMOUNT BY WHICH 
C Q CHANGES fO--’ EACH SUCCESSIVE INTEGRAnON, INTEGRATION 
C IS CARRIED OUT IN SUBROUTINE D0G12. 

C 

i)X = UX+OELX 

1F( 7.LT . 10. ) GO TO 333 
Q = 3. , ’i‘T*(Z/]3. )M>-.5.-nX 
GO TO 334 

333 g = 3.'.=»=T>;'(Z/4. -OX 

334 IF(O.Lr.EI) GO TO lDv4 
N = C 

GCTOIf 5 

lb04 li- (C. LT .OIGOTOl Jb6 
E I = E ( L ) 

XL=, iK., ( j3J1*Q 
VAL 1M=. 23^0 

CALL 0:G32(.,,XL , PCT , Z A , L 1 , 0 , T , A , 7 ) 

C AL L DOG 32 ( XL , VA L I M , FC T , ZB , E 1 , Q , T , A , Z ) 

CALL nOGG2( VAL 1M,U , FC T , Z C , E I , 0 , T , A , Z ) 

WC=ZA+Z0fZC 

WOT ( L +N- 1 ) ^ ,.0T ( L + N- 1 ) + r^QvWF ( L I 
XX( L" N-1 ) = E (L ) -Q 
J=L+N-1 

TEC D ISTi: IBUl ION FOR THF. H1ST0(,RA‘-' CuR^' ESPOND I NO TO THE 
CFiMcR or- THE BEAM DISTRIBUTION IS PlMOVEU AND STORED IN 
ViUF AND XUF AS THt UNFOLDED COUNTS AND ENERGY 
RESPECT I VELY. 

1 F ( L ,NE . ICENT I GOTQl uJ5 
WUF ( 0) -'.iD 
XljF(N)=EI-0 
K = N 

GOTO 1 . ■ 

Iv 36 CONTINUE 

no 1717 1=1,7 

ELIMINATE LOW NU*>'i3FRi;D "BINS" AS THEY ARE ONLY PAxTLY RIL- 
LED AND ARF N-JT R E PR E. SE NT AT I VE OF THE CUMULATIVE 
D I STR IBUT ICN. 

'aUE( 1 ) =WUF( 3 ) 

17 17 WGT ( I ) = ,'?T( 3 ) 

00 6 1=1,4 
6 T 3 ( 1 ) = T 6 ( n . 

GOTO 779 
777 DC77B N=1,K 



PROCESS THE UNFOLDED THEORY CURVE. SETTING EXX(l) E^OAL 
TO ZERO ALLOWS PLOTTING OR THIS THEORY CURVE WIT’-^OUT ANY 

experinental data. 



WQT (N> =WUF( r.) 

XX(N)=XUF(N) 

KS IG=I 
J = R 

775 CONTINUE 
r.XX( 1 )=D 
DO 6 1=1,4 
5 T 3 ( 1) = T 5 ( I ) 

7 79 w:orMx = . 

C01EJ7 N=3,J 

C DETER'MNE IHE MAXIMUM POINT jF THE D 1 S T R I ?U 1 J 0 . , T ' , 
C AfjD thus establish thf MOST PRUBAPLE LNE-<.GY LOSS, vM. 
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1F( HQT( K ) . LE.WQT (N- 1 ) . AND. WQT (M-1 ) .GE. WQT (N-2) ) 
IGOTOUDa 
GOTO li'wT 

1C 18 II-(iVQT(\-l) .I.T.WCTMXIGOTG 1107 
WOT MX = WOT ( N-1 ) 
f: 1 =TEMP 
'.;P=b l-XX(N-l ) 

1C- 07 CLMTIMUE 

normalize all points of niSTRlBUTION AND STORE IN HW. 

DO M=1,J 

HW(H)=WQT(M)/W0TMX 
UC9 CCNTINUE 
J = J-1 
HR I 



DLTERMINL HALFNIDIH OF THE THEORETICAL DISTRIBUTION, HVJA. 
Du L, !■ M = 1,J 

I F ( HWn + 1 ) .GT. , 3. AND.Hw(M) . L T . . 5 ) HW I = XX ( M ) + ( . 5- 
IHWIM) ) rDELX/(HW(M+l )-HW(M) ) 

1 F ( H. L-0 . J . AND. HW 1 . F Q. u ) GUTOl 1 9 
GCT0121 

119 WR ITE { 6 , ) 

hWA = D 
GoroioiD 

?u- FORMAT (/ ,?X, ' NO HALFMDTH OBTAINED') 

121 I F ( I to ( M ) . GT . . 3 . A Vi) . HW ( '• + 1 ) . L T . . 5 ) H WA = X X ( H ) + ( H W ( M ) 

1 - . 3 ) ■- Dr L X/ ( h.w ( M ) -H W ( •■’ + 1 ) ) - HW I 
U 1C CONTINUE 

l-stablish cutoff signal for plotting subroutine and 

ESTABLISH CORRECT GRAPH TITLES FOR ABSORBER USED. 

V / f j + 1 1 1 

I F ( Z . G t . A . 1 ) GO T 0 1 3 

no 2 1 = 1 , A 

2 r2(I)=Tll(I) 

GO roioi 1 

13 I F( Z.GT . 13, 1 )G 1T029 
DO 3 1^1, A 

3 T 2 ( I I = T 3 ( I ) 

GOTO] <jl 1 

29 IF(Z.Gr.29.1 )G'JT05u 

DO ! :■ I = 1 , A 
n T2( 1 )=T 1D( I ) 

Gc r 0 lull 

f'G IF(Z.GT.Sj.l)GJT06A 
DO B I =1 , A 

8 T2( 1 ) =T3 ( I ) 

Goroioi 1 

6A IF( Z.GT.tA, 1 )GuTO>32 
DC 9 1 = 1, A 

9 T2( I )=T9 ( I ) 

GCTOlOl 1 

82 DO 7 I = 1 , A 

7 T2( I ) =T7 ( 1 ) 

SEND DATA TO PLOTTING SUBROUTINE FOR PRuCESSlNG AND 

PRINT Results of ciamrutat ions in tabular form. 

1 . 11 CALL GR APH( X X , WQT , E X X , E X Y , TA , T 3 , T 2 , T 1 ,H W A , E I , T , 
lCiJTR,CTSNUR,DP ) 

WRITE (6, 22) yP,HWA,T,?,EI 
IFIKSIG.EQ.P) goto 777 
11 FORmaT(7F1 >. ■■) 

22 E'lRVAT ( / ,2X , ' : ',F7.A,' HALF WIDTH: ',F6.3, 

1' T( ''-'VS 1-CV ) : ',F6.A,' Z: 'tFA...,’ INCIDENT c.NFRGY: ’ 

2 , F . 2 , / ) 

ST IP 
END 
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SLQROUT INE D JG 3 2 ( X L , XU , FCT , Y , E I , 0 , T , 0 , Z ) 



PURPUS5; THIS SU3ROUTINE PEREORMS THE INTEGRATION REQUIRED 
IN THE P-W THEORY. 

IMPLICIT REAL*P ( A-H,n-Z ) 

A=. 500 > (XU+XL) 

H=XU-XL 

C=.4086 3 193 •92.+ 74u7R00=Np 

Y = .33(,i‘^3E5 ■> ‘^73 5..483D-24 

1 ( rci ( ( A + C) , / , D, T , E ! , Q) «-FCT ( ( A-C ),Z,D,T,EI,0)) 

C=. 49 2 ‘ns 73 3 772 A I 7D'.';=P> 

Y=Y+. 81 37197365432335D-2V 
1 ( FCT ( ( /'+0 , Z , D, T , 1 1 , 9 ) +PCT{ ( A-C ) , Z , D , T , F 1 , Q ) ) 

C=. 46238 11277937332?Dv 
Y=Y + . 1 2696 . 3263463U' 3'. 0-1* 

■1 (FCT( ( A + C) , 7 ,0, T, ri ,.l) + FCT ( (A-C ) , Z , D , T , F I , Q ) ) 

C 4674 )3J3 796 3 3')9840, -.|3 

Y = Y + . 17 1 3 6 9 i 1 4 6 6 9 1 3 7 1 7D- 1 * 

1 ( FCT ( ( A-i-C) , 7. , U ,T , E I , 'M '-FCT ( { A-C ),Z,D,T,FI,0)) 

C = . 443 1 6J5 77bB 3.,'26-.. 6U. -3 

Y = Y<-. 21 A 1 794 11113 VU 9-1* 

1 (FCT ( ( A+C) , Z , D, T , E I , QH-f CT ( ( A-C ) ,Z, U, T, E I , Q) > 

C = . 4247 830 636628499,) 

Y = Y^-.23499 '2 1631 183v 83D-1* 

1 (EC1 ( ( A. + C) , / ,9 , T,E I ,01 +FCT ( (A-C) , Z , D , T , E I , Q ) ) 

C=o 3 1724 18 9793 39 71 2 - iu *3 
Y= Y + . 2^ 3 4 2 . ■ 4 6 7 3 9 2 6 7 7 7 4 O- 1 * 

1 (ECT ( ( \ + C) , 7 ,D, T , E I , ■:) +FCT ( ( A-C ) , 7 , D , T , E I , Q ) ) 

C=. 3 66'. 9 1 JS9 Ji7 I 44 3401.: *3 

Y = Y + . 329 1 I 1 ] 13 8 6 i P 6 42 3U- I* 

1 (FCT ( ( A + C ) , Z , D , r , L- 1 , ) +ECT ( ( A-C ) , 7 , 0 , T , E 1 , C ) ) 

C = - 3 3 i. 9 2 2 1 9 3 4 6 0 1 7 6 ^ B 

Yr V4- . 1 72 WQ7,' 6 442 A 9 6 4 0- I * 

1 ( F C I ( ( 3 ^ C ) , •' , (J t T , E 1 , .1 ) + FC T ( ( A-C ) , Z , 0 , T , E 1 , 0 ) » 

C=.293H6787962..o8l 16 ) *-3 
Y=Y + , 3 9.) )3 636 163D-1* 

1 (FCT { ( A 4C ) , Z ,1) , T ,t 1 ,0) i-f-CT ( ( A-C ),Z,D,T,FI,Q)) 

C = • 2 6 3 - 4 9 4 6 4 4 6 6 1 1 4 7 i ) l; * B 

Y = Yt-. 416666621 1 -■.4733780-1* 

1 (FCT ( (,•.■^C) , Z ,0 , T ,E ( , J) <-FLT ( ( A-C ) , Z , D , T , E I , 9 ) ) 

C= . 2 K‘6 766 ^4 ,6 6 3 1767-)-. *0 
v=Yt-. 4^826 46 :j 2 1 9 . 6 i9- 1 * 

1 ( FCT ( ( 3 rC ) , 7 , U , T , E 1 , 9 ) t FC T ( ( A-C ) , Z , 0 , T , F 1 , Q ) ) 

C=. 1689 343 11 4I.-633? 

Y=Y + . 4668o rn3473Ml 9470-1 * 

1 (FCT ( ( ) , Z , T, E I , 3) +FCT ( ( A-C ) , Z , D , T , F 1 , Q ) ) 

C=. 1 1969 25 31126 .oB64 ); - o 
Y=Y-t . 46 92? 1996‘^v'4 J?28 30-l 'K 
1 ( FCT ( ( A fC ) , 7 , .9 , T , C- I , j ) +FCT ( ( A-C ) , Z , D , T , E 1 , Q ) ) 

C=. 722369a 791392260-146 

Y = Y + . 4 7 8 1 9 5 6 ,. J 3 97j3 7 4 3<- 0.- 1 * 

1 ( FCT ( ( no , 7 ,0, T , E I , 1) fFC r ( ( A-C ) , Z, U, T, E I , 0 ) ) 

C = . 2 4 1 6 3 d 32 24 3 8 6 916 go- 1 4 .3 

Y = ;3*( Yf.4o27.' I442573639r ID-l* 

1 (FCT ( ( A + C) ,7 ,L), T , E I ,0) ■^FC^ ( (A-C ) , Z , D , T , E I , Q ) ) ) 

RE I URN 
END 
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REAL FUiMCTlGN W I B ( X , E 1 , 1 , Z , A , 0 ) 



PURPOSE: This PUnCTIOM CALCULATES THAT PART OF THE ENERGY 
OISTRIBUTICN OF AN ELECTRON HUE TO LOSS OF ENERGY BY 
IONIZATION while PASSING THROUGH AN ABSORBER. 

IMPLICIT REAL--:=8( A-H,0-Z ) 

IL ARE the AVERAoE IONIZATION POTENTIALS PER ELECTRON F CR 
LEAD FUR EACH SHELL, BCGIiNNiNG V'JITH THE K SHELL. 

NPCTL IS THE NUMBER OF IONIZATION POTENTIALS USFD FOR LEAD 

IK THIS CONTEXT, L = LEAD, G = GADOL I N I U>1, 5 = TIN, C=CGPPER, 
A=ALUMINUN, AND B=BERYLL1UM. 

RfAl.NB IL('J)/. .830wU, . '1A7 j^ , 3200,. K'CAGO, 

1 .0* 7 f) , . Jv' 2 / 

• RF\L^ B NL ( G ) /2 . , 3. , n . , 32. , 18 . , A. / 

Rl-AL^-B IG(B)/.v, 20,. < -TbB , .0^1 BV , . , J02A , .. wC026/ 

H E M.* B NG (•5)/2.,B.,lB.,2B.,q,/ 

REAL^f'B IS(")/. 292,. . AA , . 0<. 72 , . 03 v 26? , . Jt COC 1 / 

R E AL-'' 8 N S { B ) /2 . , 3 . , 1 R . , 1 3 . , 4 . / 

REAL’^B IC( 3 )/. . ' 39B , .0 .) .99, . A;OC60/ 

REALM'S NC( 3 )/2 . , , 18./ 

REAL^'^a lAI 2 )/.■' ..IBt ,.2.u.uBB/ 

RTAL^^B NA(2)/2.,8./ 

KEAL-B .13(2)/, 1 h 8 , . L .,’0 1 A / 

REAL* 3 NL3(2)/2.,2./ 

NPuI L=h 
NPOTG=B 
NPUT S = E 
NP0TC=3 
NPOT A = 2 
^|;^.1y8 = p 

SB IS THE SU7'MATI0N OF IONIZATION POTENTIALS FOR THE 
EXPER INFNTAL ATOM. 

SB=\). 

XYASE IS THE REST NASS OF AN ELCCIRON IN NEV. 

XMASE=. 61 1 ->..6 

BETA IS THE MORNAl v/C USED IN RELATIVISTIC CALCULATIONS. 

GAt NA IS THE NOROAL TERM USFD IN RFLATIVISTIC CALCULATIONS 

BETA = DSQRT ( ( F I - E I + 2 . L I ^ XM A S E ) / ( E E I +2 . E I -XOASE 
1 hX -lASE^ XMASr ) ) 

GA NMA=1 . /( I .-BFTA^. -2 ( 1 ./2. ) 

at this point, the program ,MUST CU'-’PUTl- THE PROPER CON- 
ST ANTS PO:R THE EXPFR I ^'FNTAI Z. A SERIES OF LLGIC STATE- 
MENTS SELECT THE CORRECT FORMULAE. 

SEE ALUMINUM AS EXAMPLE. 

1 F ( Z.GT .4. 1 )GUT013 
DU 7 . 1=1, 'JPOTB 

•3=IB( I ) *imB( I )-OLOG( 2.v'El/( IB( I ) >i'- ( l.-BETA--^2) ) ) 

7') S3 = SB-t-i5 
B2=l>’.c4 
A 1 = . 0 o 8 I 

DELT = 4. 6C6=:-DL0G1 . { 6 E T A^^GAMMA ) -2 . 8 3 

DELT=4,6<'6X + C FCR BERYLLIUM. SEE STERNHEIXFR. 

CO TP 111 

13 I^( Z .G^ .13. ] )G0T029 

C CCMPUT2 SB for F XP E •'’ I t J T AL 2. B IS THE I0NIZ.AT1C. ’ T' - 
C TIAL FROM ONE SHELL. 



60 



I 



OOOO OOOO OOOOOOOOO 



DO 71 l=l,NPOTA 

B = IA( I )*NA( I )*nLUG( 2.>:'EI/( IA( I )«( l.-0GTA“*2) ) ) 

71 SB=S3+B 
C 

C E;STA9LISH VALUES FOR B2 AND Al. THESE ARE MATERIAL DE- 
C PENDENT CONSTANTS USED TO COMPUTE THE AVERAGE ENERGY 
C LOSS BY IONIZATION AND ARE SUPPLIED BY THE STERNHLIMER 
C PAPER. THESE ARE 3 AND A, RE SPECT I VEL Y , FROM HIS PAPER. 
C 

82=16.77 
Al=.C7AO 

DELT IS A PART OF THE EXPRESSION BY STERNHEIMFR FOR THE 
AVERAGE ENERGY LOSS DUE TO IONIZATION. NUMBERS HEREIN 
REPRESENT OTHER STERNHEIMER CONSTANTS AS FOLLOWS: 

SMALL A .0905 
SMALL M 3.5U 
XI 3. LUO 

-C 4.210 



DELT =4. 606YDLOG10( BE T A«G AMM A ) -4 . 2 1 + . C 906* 

1 (3.-DL0G1C(BETA*GAMMA) )**3.51 
GO TO 111 

29 IFIZ.GT.29. 1)GOT050 
DO 72 I = 1,NP0TC 

B= IC( I )*NC( 1 )*DLQG(2.*EI/{ IC( I )*( 1.-8ETA**2 ) ) ) 

72 SB=SB+B 
62=15.09 
A1=.07C1 

DELT=4. 5C6*DL0Gin( BE TA*GAMMA )-4. 74+. 119C* 

1 (3.-DL0G10(BETA*GAMMA) )**3.^8 
GO 10 111 

50 IF( Z.GT .50. 1 )G0T064 
DO 73 I=1,NP0TS 

B = IS(I)*NS(I )*DLOG( 2.*EI/ ( IS( I )*( 1. -BET A** 2) ) ) 

73 SB=SB+B 
B 2= 13. 83 
A 1 = . 0 64 7 

DELT=4. 606*DLOGIO( B E TA*GAMM A ) -6 . 28 + . 404* ( 3. - 
IDLOGIOI BGTA*GAMMA) ) 

1**2. 52 
GO TO 111 

64 IF(Z.GT.64.1)GQT082 
DO 74 I=1,NPGTG 

B= IG( I )*NG( I )*DLOG( 2.*EI/( IG( I )*( i.-BETA**2 ) ) ) 

74 SB=S3+B 
B2=13.3 
Al=.0624 

DELT=4.606*DLOGlj(8ETA*GAMMA)-6.8+.418*(3.- 
IDLGGIOI BE TA*GAHMA) ) 

1 * * 2.10 
GG TO 111 

82 IF(Z.GT.82.1)GGT0112 
DO 75 I=1,NPQTL 

R=IL( I )*NL( n*DL0G(2.*EI/( IL( I)*( l.-BETA**2) ) ) 

75 S3=SB+B 
62=12.81 
Al=.0608 

DELT=4. 6C6=^DL0G10( BETA*GAMMA ) -6 . 9 3 + . 0 652* 

1 (4.-DL0G10( RETA*GAMMA ) ) **3.41 
GO TO 111 

AR IS THE QUANTITY SMALL A, FROM BLUNCK AND WESTPHAL, 
MULTIPLIED BY R, THE TARGET THICKNESS IN CENTIMETERS. 

Ill AR=U. 154*Z*T/A/BETA/BETA 

CCMPUTE BS2, THE SMALL B SQUARED FROM THE BLUNCK AMD 
LEISEGANG PAPER. 

BS2=3.*S3/ ( Z*AR) 
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C VAKl, VAP?, AND VAR? ARE COrjVEr,’ I tAiC E STORAGE LOCATIONS 
C UStn TG STORE CP^^PUTTr) PORTIONS OF S rERNHE 1 MER ' S AVERAGE 
C ENERGY LOSS. 

C 

VAR1 = A1-T/BTTA>.'*;? 

VAR2 = B2 + .4? ♦■Z. ’"DLOG ( BETA<--GAMNA ) 

VAR3=DlCG(EI )-3ETA^-2 

QAVE IS THE FINAL VALUE FOR THF AVERAGE ENERGY LOSS. 

QAVF=VAR 1-M VAR2 + V AR3-DELT ) 

B1 =US0RT ( BS? ) 

X IS THE PORTION OF SOKE TOTAL ENERGY LOSS Q DUE TO 
IONIZATION. 

X = Q-X 

XLAMB IS THE LAKBOA VARIABLE FROM THE BLUNCK AND WESTPHAL 
P A P E R . 

XLAMu= ( X-QAVl ) / ARi-DLOG( EI/AR )-l . 116 
WUJM=0. 

TEE variable TEST IS USED TO KESTi’ICT OUTPUT TC NUMBERS 
OF LSEEUL MAGNITUDE. 

TEST = Xl ( B 1 B U ? . 24 ) 

I '■ ( TE S T . GT .133.) GO T 0 1 1 2 

CCEEF ARE THE CE E F F I C 1 E NT S FOR EACH OF THE FOUR GAi;SSIAN 
CLI-’VES BLUNCK AN') NESTPHAL USE TO FIT I HE LANDAU UISTRI- 
3UT10M. AN ADOn lONAL CONSTANT HAS BEEN MULTIPLIED INTO 
EACH FOR CONVt.O ENCE IN RUNNING THE PROGRAM. 

rnt-PF 1 = ( M , ] V 1 . B ) / OSCk I ( B ] - B ) + 1 . a -U . b ) ^ , iU / 

COEFE2= ( . ) /ns ART ( B 1 4-2.E *2. ) *. 1 D7 

CHEEP 3 = ( j . • 1 ■)'f 3 . ) / ■'SDR T ( B 1“!^ B 1 + 3. 0‘'^3 , P . 1 D7 

CUE ="04 = ( . .■• ' .) / USQRT ( B 1 v B I + s . o . IL'7 

GAOS ARE THE B-A GAUSSIAN CUP.VES INCLUDING C OF. E F I C 1 E N T S . 

GAUSl=COPFEP'0'fXP( -(XLAMB-' . . ) *n« 2/ ( B1 *3 ) ♦•I .3=^1.H) I 

G A'J S 2 = C 0 c F F 2 " D E X P ( - ( X L A B - 3 . l M 2 / ( B P- B 1 + 2 . J 2 . ) 

GA'JS3 = C JEFF3*DEXP( -(< L AM B-6 . B ) "•>: ?/ { BP^B 1 • 3 . >---3 . ■ ) ) 

G \USA = C0EFF4>:-Dl. XP( - ( xL \MB~1 1. .) ) -M-Z/ { Bl PP ) ) 

CORINT IS AN EMPERICALL^' DERIVE'.) CORPECTIO.s FACTOR T'O C.Pi- 
PEnSATE rOR TF'E !-A(,T 1 HA I THE 0 1 S I P PUJE 1 NV FUNCT 1 uN C-N- 
NOl BE INTEGRATED EXACTiY AS THE RAniAEION LESS 1 S I ' 1 - 
BUriON BFCCMFS INFINITE AS ITS VARIABLE (ENERGY LOSS) 

approaches zero. 

CORINT=l .2CM ( 1 i ,32 3 . ^ 3 ) 

GAUSl = GAUSP<Cu U NT 

V;iON IS THE SUN IF GAUSSIAN CURVES REPRESENTING THE ^^1ERGY 
DISTRIBUTION DUE TO lOOiZATlL-N LOSSES. 

W I (J N = G A U S 1 + G A U S 2 + G A U S 3 + G A U S 4 
112 RETURN 
END 
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RtAL FUNCTION WP AD 8 ( X , f. 1 , T , Z , A , 0 ) 

PUPPUSS: THIS FURCTIJ",! CALCULAF'^S THAT PART CF THE ENERGY 
DI STR ir.UT ! CN OF AN FLECTROivj DUE T C LOSS OF ENERGY BY 
RACIATION WHILE PASSING THROUGH AM ABSORBER. 

IMPLICIT REAL*8 ( A-H,')-Z ) 

U IS AN INTL-RMEUI ATE STORAGE LOCATION USED IN COMPUTING 
ALPHA. 

U=183.v-/Z*«{ 1. /B. J ) 

CCf'PUTE VALUE OF ALPHAxR FOR TARGET THICKNESS, ATOMIC 

NUMBcR AMD ATOMIC WEIGHT OF EXPFR I v,Fi>^T AL ABSCRBEP. R IS 
THE TARGET THICKNESS EXPRESSED IN CENTlMcytRS. 

ALPHP = . . j01A'NT!=ZvZ/ A M 4 . .■/ 3 . v *DL0G ( U ) + 1 . /G.v ) 

CCPPLiTE VALUE FOR B, A NMPMaLIZIrG CONSTANT. 

B^1./1)GA.MMA( ALPHK+1 . ) 

CCiMPUTE „RAD, THE NUMBER OF COUNTS EXPECTED AS A RESULT OF 
PACIATin,'; energy LTSSE'S. X RcORESRNTS THAT PORTIuN OF 
Sl*-’E total ENcFGY LjSS WHICH IS ..01 LOST B'' 1 O.n, I Z A T I UN . 

WRAD= ( B*ALPHR>^ ( X-=TA LPHR I / ( E I M^^ALPHR ) ) / X 

RETURN 

END 



f? F AL FU iNCT I ON BE AM* 8 ( E 1 , HW , X ) 

PURPOSE: THIS FU.NCTIDN REPLICATES THE M S T R I ''UT I ON OF 
'iEAI ENE.-vGY as a SIMr’LE GAUSSIAN CURVE BASED ON THE 
EXPERIMF Nil ALl Y OBSf RVE J HALEVIDTH AN'^ I S LGCAIED 
RcLATIVf: TO THE BE AH E.N'ERGY cjY THE CJMPUTlO AVERAGE 
RECOIL LOSS. 

i:'PLICI T RE AL*8( A-H,G-Z) 

Y = 2. 

ALFA = 4,*DL0G( Y )/( Hi-.*HW) 

I3cA.V = Df XP(-( (X-E I 1**2) J^ALFA) 

RPTURN 

END 



REAL FUNCTION FC T* R { X , Z , A , T , E I , 0 ) 

PURPOSE: THIS FU.ACTICN ^'ULT 1 PL I CAT 1 VF LY JOINS TOGETHER 1 HE 
PREDICT 1L-! OF ENERGY' LOSS DUE TO IONIZATION AND Th£ LOSS 
DUE TO RADIATIO.N. 

IMPL ICI T RE AL*B( A-H.O-Z ) 

FCT = '-;R \0 ( X, E 1 , T , Z, , 0 ) '•■.!’ N. ( , E 1 , T , Z , A , 0 ) 

RETURN 

E.\0 
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SUB ROUT I ME GR A Pri ( TH X , T H Y , EXX,EXY,T1,T2,T3,T4,HW,E1,T 
1CURR,CTSMCR,QP ) 

PURPOSE: THIS SURRUUTINF PLOTS THE RESULTS OF THEORY CAL 
CULATIGNS AND FXPLRIPEMTAL RESULTS. 

I'-IPLICI T REAL^"8( A-H,0-Z ) 

REAL>J=d TITLEX(l) 

DATA TITLEX/'MEV '/ 

REALMS T ITLL Y( 1 ) 

DATA T I TLEY/‘Ni1RMAL IZED COUNTS '/ 

RFAL<=8 TITLEUI?) 

cata t I tleu/ ■opp'^dahl box 0 •/ 

DIME MS 1 ON T1 (4 ) , T2{ , T3( A ) , T4( A ) , THX ( 3CC ) , THY( 3CC ) , 

iL XX(2U. ) ,EXY(?'. ' ) , Cl SUP ( 2UL ) ,CI SDN( 20 ., ) 

DETEHMUiE SCALING FACTORS A\D CHANGE THEORETICAL CURVE 
TO PR0'E;R LUCA! I UN. 

X = 8. 

Y = S. 

S1XIN=( E I-0P)^-2, 

IS1XIM=S1X1N+. 3 
SI XIM2=I SI XIM 
SI XiN=S I XIM2/2. 

If (HN.I-Q.L’)GOTT<) 

X^.CR2=; 

INCR? = XNCR24 
IF{ 1MCR2.lt. 1 ) INCK2 = 1 
X I NCR 2 = 1 NCk 2 
XINCK = XlMC‘^2/2. 

SP=S 1 XI N-tj, -i-XI ,\CR 
It-(SP.LT.O)GOTO'^.. 

GO TOPI 
SP = '' 

INCR=( EI-1.3AT )/6.4.5 
XINC‘^ = INCR 
91 N = C 
1 = 2 

20 1=1+] 

]F(LXX( 1 ).GT.l. ) G0TC21 
GO TO 2 2 

21 N=N+] 

GOTO 20 

22 IF(lj., E0.'j)i;0T0 1.; 

DO 23 J=1,M 

E XX ( J ) = ( E X X ( J ) - S P ) / X I NC R + C OR'? / X 1 NCR 
IM J.EO.N) GCTf}23 - 

23 CCNTINUE 
1 ;C 1 = 3 

= - ) 

30 I = i + 1 

IF { THX( I ) .G1 . 1 . ) GOTO 31 
GOTU22 

31 M=M+1 
G0T030 

22 DO 33 J=1,M 

THX( J) =( THX ( J ) -Si^) /XI NCR 
IF( J.b0.M)GCr033 

IF( THY( J+1 ) .GT .THY ( J) ) I HYMAX= THY( J+ 1 ) 

23 CCNTINUE 
IF(N.EQ.0)GCT01C I 
SCALE=A. 

IFIT.GT .A. 3) SCALE=3. 
no A) J=1,N 

LXY( J) = EXY( J )/CTSMtjR"^SCALE 
AG CfilTINUE 
1 - 1 00 5 J = l, 

SCALE=a. 

IF(T.GT.A.P) SCALE=3. 
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THY(J)=THY(J) / THYMAX'-'SCAL t 
50 CONTINUE 
C 

C SORT EXX AND ELIMINATE VALUES OFF GRAPH. 
C 

IF(N.EQ.U) GOTO 102 
K = N-1 

DO 60 1 =1 ,K 
IP1=1+1 

00 6U J=IP1,N 

IF(tXX( I ) .LE.EXXU ) ) GOTO 6u 
TEMP-'FXX ( I ) 

TEV=EXY(I) 

EXX( I ) =EXX( J> 

FXY( 1 ) =f XY( J ) 

E X X ( vl ) = T E ’-1 P 
EXY{ J ) =TEM 
5) CONTINUE 
J= .) 

K K = 0 
K = N 

on 6 1 . i = i,;< 

1F(KK.E0.1)GUT061 
IF( EXX( I ) .LT.G )GUT06l 
J = J<-I 
N = J 

EXX( J ) =EXX( I ) 

E^YI J) =EXY( i ) 

1 F ( EXX( J ) ,GT. X ) GUT062 
GOTO 61 

62 K=J-1 
KK^ I 

61 CONTINUE 
IU2 KK=0 

J=v; 

»< =H 

DO 71 l=l,K 
IFIKK.LQ. l)GCm71 
IF(ThX( I ) .LT.'j )G0Tn71 
J = J+1 
M = J 

TFX( J ) = THX ( I ) 

ThY( J ) = TH.Y( 1 ) 

IF( THX{ I ) .01 . X )GUT 172 
Gcrn^i 
72 M=J-1 
KK=l 

71 CONTINUE 
C 

C DETERMINE ERROR GA MAG.NITUDE. 

C 

IFIM.EQ.t.;) GClOl j3 
DO 7 I=1,N 

CTSOMI I )=EXY( I )-lSURT ( rXY ( I )*2.-. . ) / 2 
CTSUP( I ) =EXY( I ) +USORT lEXY ( I ) ”2 j;'. ) /2 
7^,’0 CONTINUE 
C 

C INITIALIZE PLOT AND XRITE IDENTIFICATION 
C 

1>.3 CALL PI ITS 

CALL SYMIOLI , , . 2 R , T I TL EU , u , 1 6 ) 

CALL PLOT r , 2., -3) 

C 

C DRAW CJUTLir.E OF GRAPH. 

C 

Z = G' 

CALL PLnT(7,Y,2) 

CALI RL0T(X,Y,2) 

CALL PL0T(X,/,2) 

C'.LL PLUMZ,Z,2) 

C 

C DRAW OUTLINE OF TITLE BGX. 
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Bl = . 

B2=4.4 
B3 = 3. 

B^* = 4. 6 
B5 = B1 + . 1 
B6 = B3+-. 1 

CALL PLOKBl ,R-^ , 3) 

C ALL PL )T ( B2, R3, 2) 

CALL PL0T(B2,BA, 2) 

CALL PLOT(« ] ,BA,2) 

CALL PLH r ( B1 , R3, 2 ) 

CALL SYP30L ( -35 » Bo, . lA , T1,0,C,32) 

B5 = P.5 + . 24 

CAl.L N.JMBEK(B5 , D6, . 14,T, },2,3) 
B5=B5+1 . 32 
Bo=BA+» 1 
EXP=2. 

CALL f-jUMdER (B5 , Bo , , 37 ,P XP , V* V. , -1 ) 
f»5=B5 + l , 20 
Bb = BC'— g 1 

call UU":dER ( B5 , Bo, . 14 , P I ,0,0,2) 
B5=B5-2 . 76 
B6=B6+. 24 

CALL symbol ( B5 , 36, . 14, 

I p ( rj . b 0 . .. ) GO ro 7 I 
36=B6+.J7 . ■ 

R5-B6+. B2 ■ • 

B7=B5+. 5d 

CALL PLOT (B5,Bo,3) 

CALL PLOT (B7,B6,2) 

B 6 = B 6~ • jl 
B5=B3-. 82 
7)1 Bo=B6+,24 

CALL SY'U50L (33 , 36, . 14, 

R/S = R.<',<- . R A 

CALL S Y M HQ L ( '’ 3 ; B , 1 4 , 

PLOT TIC -MRKS. 

G=. 1 
H — V” 1 

15 7'j CALL*PL0T(G,Y,3) 

CALL PLOT ( 0,1 1, 2) 

G = G + *. 1 

IF (G.LE.X)GO TO 137o 
G = X-. 1 
1 

158v, CALL PL0T(G,M,3) 

CALL PLOT (G , . , 2 ) 

G = G~ • 1 

1 Fi’gIgT.O. )GU TO 13 3Q 
G = 1 . 

H “ 2 

167o CALL PI OT (G,Q, 3 ) 

CALL PLOT (G , 2 ) 

G=3+l. 

IF(G.LT.X) GO TO 1670 
H=Y-,2 

1680 CALL PLOT(G,H, 3 ) 

CALL PLJT(G,Y,2) 

G~G~ 1 < 

fPl’G.GT.' . ) GO TO 1 63 j 
G=. 1 
H — Y"“ 1 

1371 CALL* PLOT (G,H, 3 ) 

CALL PLOT ( . ,M ,2 ) 

H=H-. 1 

IFIH.GT. j. ) GO TO 1 57 1 
G=X-. 1 
1^=. 1 

1581 CALL PIGT(G,M,^) 

CALL PL0T(X,H,2) 



TZjO,-^, 32) 



T3 ,0.) , 32 ) 
T 4 , 0 . i ,321 
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H = H + . I 

n-(H.LT.Y) GU TU 15B1 
G=X-.2 

1671 CALL PL0rio,H,3) 

CALL PLUT(X,H,2) 

H=H-1 . 

IF(H.GT.J) G!) TO 1671 

G=.2 

H=l. 

1681 CALL PinT(G,H,3) 

CALL PLOT( ,H, 2 ) 

H=H+1 . 

IF(H.LT.Y) GO TO 1681 

label axis. 

CALL SYMBOL! ^ .6,-. 6 lA ,TI TLEX,C 

CALL SYMBUL(-. 26 , i .56, . lA, Tni.EY,9.i. , 24) 

NUMBFH EMEGGY AXIS. 

1.5-' G = -.2B 
H=-. 32 
FLO=SP 

lc6.. CALL iM'JMBFR(G,H, .14,FLU,^',2) 

G=G+1. 

IF(G.Gr.X) GO TO 1.7 
FLJ = F HH-Xi 'vC.\ 

GO TO lv<60 

PLOT EXPEK O'lENT AL DATA. 

167J IF(N.EO.J)GOrOIu4 

C A I I, L I OF ( F X X . F X Y , ,0 .2,-6) 

PLOT ERROR FLAGS, 

DO 11 1=1,0 

ESl = rXX( I )-o 5 

ES2 = EXX ( 1 ) -I . 5 

CALL PLOT(LSl,fTSUP( I ) ,3) 

call plot ( c.S. , CTSUP ( I ) , 2 ) 

CALL PLOT (i:XX( 1 ) ,CTSJ^"( I ) ,3 ) 

CALL PLOT (r. XXI I ) ,C I S0M{ I ) , 2 ) 

CALL PLOT ( F S] , CTSD m ( I ) , 3 I 
call pi OT( f S2,CTSD0( I) ,? ) 

1 It; CCOriNUE 

PLCT THEORY CURVE, 

104 CALL LIME( FhX, THY,M,2, 1 ) 

FMALIZE PLOT. 

CALL PLOT! .,8., -3) 

CALL PLOTE 

RETURM 

END 
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